Polaron physics and crossover transition in magnetite probed by pressure-dependent 

infrared spectroscopy 
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The optical properties of magnetite at room temperature were studied by infrared reflectivity 
measurements as a function of pressure up to 8 GPa. The optical conductivity spectrum consists of 
a Drude term, two sharp phonon modes, a far- infrared band at around 600 cm -1 , and a pronounced 
mid-infrared absorption band. With increasing pressure both absorption bands shift to lower fre- 
quencies and the phonon modes harden in a linear fashion. Based on the shape of the MIR band, the 
temperature dependence of the dc transport data, and the occurrence of the far-infrared band in the 
optical conductivity spectrum the polaronic coupling strength in magnetite at room temperature 
should be classified as intermediate. For the lower-energy phonon mode an abrupt increase of the 
linear pressure coefficient occurs at around 6 GPa, which could be attributed to minor alterations 
of the charge distribution among the different Fe sites. 
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I. INTRODUCTION 

Magnetite (Fe304) belongs to the major class of highly 
correlated electron systems and shows many interesting 
physical properties, which triggered a large number of 
experimental and theoretical investigations over the last 
60 years. Fe3 04 has an inverse cubic spinel structure at 
ambient conditions and is in a mixed-valence state de- 
scribed as [Fe 3+ ]A[Fe 2+ +Fe 3+ ] B 04, where A and B de- 
note the tetrahedral and octahedral sites, respectively, 
in the spinel structure AB2O4, with space group FdSm. 1 
Most of the proposed conduction mechanisms are based 
on either band or polaron hopping motion of the ex- 
tra electron through the octahedral sites occupied by 
Fe 2+ and Fe 3+ ions^— Magnetite undergoes a Verwey 
transition at T v «120 K at ambient pressure towards 
an insulating stated On cooling through the transition 
temperature, the dc conductivity drops by two orders of 
magnitude concurrent with a first-order structural phase 
transition from cubic to monoclinic symmetry^— The 
presence of the mixed-valent Fe ions on the B sites mo- 
tivated Verwey and others^ to hypothesize that above 
T v electronic exchange takes place whereas below T v 
the transition is related to charge ordering (CO) in the 
mixed-valent oxide. Recent resonant x-ray diffraction ex- 
periments have specified the charge disproportionation 
for T<T V as [Fe 3+ ] A [Fe 2 - 5+,5 +Fe 2 - 5 - ,5 ]B04, with S 12 = 
0.12±0.025 for one kind of Fe atom and <S 34 = 0.1±0.06 
for another kindAi~— 

On the contrary, based on Mossbauer and powder x-ray 
diffraction studies^— it was claimed that while cool- 
ing down at ambient pressure a coordination crossover 
(CC) takes place at around T v , where the spinel struc- 
ture changes from inverse to normal: 

[Fe 3+ U[F e 2 - 5+ +Fe 2 - 5+ ] B 4 -> [Fe 2+ } A [Fe 3+ +Fe 3+ } B 4 
Within this scenario, Tec increases with increasing pres- 



sure and reaches room temperature at «6 GPa, while 
T v decreases with increasing pressure ! 13 ' 14 Concomitant 
with the CC transition a small change in the pressure de- 
pendence of the tetrahedral and octahedral volumes was 
observed at around 6 GPa at room temperature, but no 
resolvable change in the spinel-type crystal structure or 
unit cell volumei 14 ' 16 ' 17 Furthermore, the hyperfine in- 
teraction parameters as a function of pressure show an 
anomaly at «7 GPa, associated with a pressure-induced 
discontinuous change of the Fe-0 bond length.!- The pro- 
posed valence transition from inverse to normal type is 
sluggish in character and hence spread over the pressure 
range 7-15 GPa at room temperature^ 

The proposal of the occurrence of a CC transition 
in magnetite at room temperature has been recently 
discussed based on a pressure-dependent thermoelectric 
power study combined with electrical resistance and sam- 
ple's contraction measurements^ and a high-pressure x- 
ray magnetic circular dichroism study* 1 ^ The results of 
both investigations exclude an inverse-to-normal spinel 
transition up to 20 GPa. Furthermore, in Reffla a 
new pressure-induced crossover near 6 GPa was pro- 
posed, with either a pressure-induced 'perfection' of the 
electronic transport via mobile polarons, or a pressure- 
induced driving of the dominant inverse spinel configu- 
ration to an 'ideal' inverse spinel. 

In the present work we carried out infrared reflectiv- 
ity measurements to study the electronic and vibrational 
properties of magnetite as a function of pressure at room 
temperature. The goal of our study was two-fold: first, 
to confirm the polaronic character of the electronic trans- 
port by monitoring the typical spectral signatures of po- 
larons in the infrared frequency range as a function of 
pressure, and second, to test the proposed scenarios for 
the crossover transition at ~6 GPa based on the changes 
in the spectral features at around this critical pressure. 
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FIG. 1: Dc resistivity p as a function of temperature. Inset: 
Arrhenius plot and fit (dotted line) of ln(Tcr) between 250 
and 300 K according to Eq. (JS| for small polaron hopping. 



II. EXPERIMENT 

Single crystals of magnetite used in this work were 
prepared from a-Fe203 using a floating-zone method 
with radiation heating ^ The quality of the crystals was 
checked by electrical transport measurements, showing 
a sharp increase of the resistivity by a factor of about 
100 at the critical temperature T v ~ 122 K, which is 
characteristic for the Verwey transition (see Fig.[T]). Fur- 
thermore, the room temperature value of p compares well 
with earlier reports i 4 ' 21 ! 22 

In the pressure-dependent studies a clamp diamond 
anvil cell (Diacell cryoDAC-Mega) equipped with type 
II A diamonds, which are suitable for infrared measure- 
ments, was used for the generation of pressures up to 
8 GPa. Finely ground Csl powder was chosen as quasi- 
hydrostatic pressure medium. The single crystal was pol- 
ished to a thickness of ss50^m and a small piece (about 
140 fim x 140 /xm) was cut for each pressure measure- 
ment and placed in the hole of the CuBe gasket. The 
pressure in the diamond anvil cell was determined by the 
ruby luminescence method. 23 

Pressure-dependent reflectance measurements were 
conducted at room temperature (RT) from 200 cm -1 to 
18000 cm -1 using a Bruker IFS 66v/S Fourier transform 
infrared (FT-IR) spectrometer. To focus the beam on 
the small sample in the pressure cell, an infrared mi- 
croscope (Bruker IRscope II) coupled to the spectrom- 
eter and equipped with a 15x magnification objective 
was used. The optical measurements were partly car- 
ried out at the infrared beamline of the synchrotron ra- 
diation source ANKA, where the same equipment is in- 
stalled. Reflectance spectra were measured at the inter- 



face between the sample and the diamond anvil. Spectra 
taken at the gasket-diamond interface served as the ref- 
erence for normalization of the sample spectra. Vari- 
ations in the synchrotron source intensity were taken 
into account by applying additional normalization pro- 
cedures. The reflectivity spectra are affected at around 
2000 cm -1 by multiphonon absorptions in the diamond 
anvils, which are not completely corrected by the nor- 
malization procedure. This part of the spectrum was 
interpolated based on the Drude-Lorentz fitting. All re- 
flectance spectra shown in this paper refer to the absolute 
reflectance at the sample-diamond interface, denoted as 
R s _d ■ The reflectance spectrum of the free-standing sam- 
ple (not shown) was found to be in good agreement with 
earlier results 4^* 

The real part of the optical conductivity, which directly 
shows the induced excitations in a material, was obtained 
via Kramers-Kronig (KK) transformation: To this end 
the measured reflectivity data were extrapolated to low 
frequencies with a Drude-Lorentz fit, while data from 
Ref|j were merged at high frequencies. Moreover, the 
sample-diamond interface was taken into account when 
performing the KK analysis, as described in Refl25l. 



III. RESULTS 

The infrared reflectance spectrum and the correspond- 
ing real part of the optical conductivity of Fe3C>4 as a 
function of pressure at RT are presented in Fig. [2j With 
increasing pressure the reflectance increases monotoni- 
cally for frequencies below 5000 cm^ 1 (w 0.6 eV), indi- 
cating a growth of spectral weight within this frequency 
range. In contrast, the reflectance is basically pressure 
independent above 12000 cm _1 (w 1.5 eV)2& We ob- 
serve two Ti u oxygen phonon modes near 355 cm -1 and 
565 cm -1 , consistent with earlier reports^ The higher- 
frequency mode is attributed to the stretching of the A- 
O bond, whereas the lower-frequency mode is related to 
the motion of the oxygen perpendicular to this bond i 27 i 28 
The real part of the optical conductivity of Fe3C>4 is pre- 
sented in Fig. [5] (b) for various pressures. The finite con- 
ductivity in the far-infrared range suggests a Drude-type 
contribution, indicating a metallic-like character of mag- 
netite already at low pressure. There is a pronounced 
mid-infrared (MIR) absorption band located at around 
5000 cm~ 1 (sa 0.6 eV), which is a typical signature of 
polaronic excitations^ - — and was attributed to the in- 
duced hopping of polaronic charge carriers between the 
Fe 2+ and Fe 3+ ions on the octahedral sites in magnetite. 4 
An additional far-infrared band is observed in the opti- 
cal conductivity spectrum at around 600 cm -1 , whose 
nature is discussed in section IIV Al 

To better quantify the changes of the optical conduc- 
tivity with increasing pressure, we fitted the conductivity 
spectra with the Drude-Lorentz model. The MIR band 
was modeled by two Lorentzian functions and the far- 
infrared band at around 600 cm -1 was described by one 
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FIG. 2: (a) Pressure-dependent reflectance at RT and (b) 
pressure-dependent real part of the optical conductivity at 
RT obtained by KK analysis. The arrow in (b) marks the 
far-infrared band as discussed in the text. The small features 
at around 3600 cm - are artifacts due to multi-phonon ab- 
sorptions in the diamond anvil. 



Lorentzian oscillator. As examples, we show in Fig. [3] 
the fitting of the optical conductivity <j\ at P = 1.5 and 
7.6 GPa including the various contributions, i.e., Drude 
term, two phonon modes, the far- infrared band, the MIR 
band, and higher-energy interband transitions. This way 
we were able to extract the contributions for each pres- 
sure applied. 

With increasing pressure the MIR absorption band 
slightly narrows and redshifts, while its spectral weight 
[SW = f o"i (uj')dui'] decreases with increasing pressure 
(see inset of Fig. [4j. The observed redshift confirms 
the polaronic nature of the band: According to polaron 
theory 3 ^ the frequency of the polaron band is a measure 
of the polaron binding energy and thus of the electron- 
phonon coupling. In general, the electron-phonon cou- 
pling tends to decrease under pressure as a result of the 
combined band broadening and stiffening of the crys- 
tal lattice. Hence, one expects a decrease of the po- 
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FIG. 3: Fit of the RT optical conductivity <n at (a) 1.5 and 
(b) 7.6 GPa using the Drude-Lorentz model. Inset: Plasma 
frequency up of the Drude term as a function of pressure 
obtained from the Drude-Lorentz fitting. 



laron binding energy under pressure. The pressure- 
induced redshift of the polaronic absorption band was 
also observed for other transition metal oxides i 30 ' 37 The 
evolution of the far-infrared band positioned at around 
600 cm -1 with pressure is depicted in Fig. [5] It shifts 
monotonically to lower energy and its spectral weight in- 
creases, with a saturation above ~6 GPa (see inset of 
Fig. [5]). 

The stiffening of the lattice under pressure, as men- 
tioned above, is demonstrated by the hardening of the 
two phonon modes (see Fig. [6]): Both modes shift linearly 
with increasing pressure. The corresponding linear pres- 
sure coefficient B was obtained by fitting the peak posi- 
tions with the function ui(P) = A + B ■ P, where P is the 
applied pressure. For the higher-energy mode we obtain 
5=4.9 cm _1 /GPa for the whole studied pressure range. 
For the lower-energy mode the linear pressure coefficient 
is B=l.9 cm _1 /GPa up to «6 GPa; at this pressure an 
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FIG. 4: MIR band as a function of pressure obtained by 
Drude-Lorentz fitting of the optical conductivity spectra. In- 
set: Spectral weight of the MIR band as a function of pressure. 



abrupt enhancement of the pressure-induced hardening 
occurs, with i?=4.9 cm _1 /GPa for the pressure range 
6-8 GPa. Interestingly, at the same pressure («6 GPa) 
an anomaly is found in the pressure dependence of the 
tetrahedral and octahedral volume,— the hyperfine in- 
teraction parameters^ and the sample's contraction. 1 ® 
In addition, drastic changes are observed in the pressure 
dependence of the thermopower and electrical resistance 
near 6 GPa.-^ This issue is discussed in section PV Bl 



IV. DISCUSSION 
A. Polaronic excitations 

In order to substantiate the qualitative picture of the 
polaronic transport in magnetite, we have analyzed our 
optical conductivity data using various theoretical po- 
laron models. Since the nature of the polaron is reflected 
by the shape of the MIR absorption bandj^i we fitted the 
MIR band in magnetite using (i) a large polaron (LP) 
model and (ii) small polaron (SP) models. 

(i) For large polarons the optical conductivity due to 
the photoionization of the charge carriers as obtained by 
Emin 3 -^ is described by 



<7i(w) = rip 



Me 2 1 {k{u)Rf 



3 m uj [1 + (k(uj)R) 2 



(2) 



where np is the polaron density and R the polaron ra- 
dius. The photoionized carrier is treated as a free par- 
ticle with mass m and wavevector k. The wave vector 
is defined as k — yj2m(Tiui — 3E p )/h, where E p is the 
polaron binding energy or the ground state energy of the 
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FIG. 5: Far-infrared band for selected pressures, obtained 
by Drude-Lorentz fitting of the optical conductivity spectra. 
Inset: Position and spectral weight of the far- infrared band 
as a function of pressure. Lines are guides to the eye. 



polaron. The threshold energy for the photoionization 
of a large polaron is 3E p £^ By assuming the effective 
mass for the large polaron to be equal to the free elec- 
tron mass m e , the best fit of our data by using Eq. ([2]) 
was achieved with the parameters np — 3.0- 10 21 cm~ 3 , 
E p = 327 cm" 1 (41 meV), and R = 1.3 A. The LP 
model describes the asymmetric lineshape of the MIR 
band well, as illustrated in Fig. [7] Assuming that the con- 
duction via polarons takes place among the octahedral B 
sites of the spinel structure, we can estimate the effec- 
tive number of charge carriers, N e ff, from np according 
to N e ff=npxV/N, where N=16 is the number of B-sitc 
atoms and V^=(8.394 A) 3 is the volume of the fee unit 
cellJaL The resulting iV e //=0.11 can be compared with 
the effective number of carriers per site obtained from 
the spectral weight of the MIR band for LPs (see Fig. [7]) 
according to A e // = [2m e V/(ire 2 )]/N J ai(uj')duj' , where 
m e is the free electron mass, and V the unit cell volume. 
The so-obtained value N e f f =0.24 is about a factor of two 
higher than the value from the fitting of the MIR band. 
Furthermore, the value of the LP radius R is smaller than 
the B-B distance, a c s/2 / 4=2 .97 A with the lattice con- 
stant a c =8.394 A ) 11 ' 35 which is in disagreement with the 
LP scenario, where the LP radius should extend over mul- 
tiple lattice sites M- This questions the applicability of the 
LP model for describing the optical data of magnetitej 3 - 

(ii) The optical properties of small polarons have re- 
cently been reinvestigated in the frame of the Holstein 
model using the dynamical mean- field theory^ In ReflU 
both the antiadiabatic and the adiabatic regimes for the 
polaron formation are discussed, which are distinguished 
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FIG. 6: Phonon frequencies as a function of pressure obtained 
by Drude-Lorentz fitting of the optical conductivity spectra. 
The lines are linear fits to the data points, in order to obtain 
the linear pressure coefficients. 



by the adiabaticity ratio "f=uJo/D, where loq is a typical 
phonon frequency and D is half the free-electron band- 
width. We can estimate from our optical data the value 
of the adiabaticity ratio to 7«0.1 by using the theoretical 
value 39 £>«0.7 eV (since we have wo~75 meV). Hence, 
we can refer to the formula in Ref.38 obtained in the 
adiabatic regime. 

Furthermore, the phonon-induced broadening of the 
electronic levels, given by the variance s of the phonon 
field, compared to the electronic dispersion represented 
by D, is crucial for the shape of the polaronic optical 
absorption^ The variance s is given by 3 - 8 - 



l (T) = E P u Q cothuJo/{2T). 



(3) 



We estimate the value of Ep from the position of the 
MIR band, which is located at ui max ~2Ep for small po- 
laron excitations ) 34 ' 38 and thus obtain i?p«366 meV. 
With the typical phonon frequency of a;o~75 meV for 
magnetite this gives s~165 meV at room temperature. 
Thus, the phonon-induced broadening is of the same or- 
der of magnitude as the electronic band dispersion with 
£>«700 meVi^S Therefore, we will consider both limits 
(s^>D, s<C-D) in the following. 

For s^D the absorption by small polarons is described 
by a Gaussian-shaped peak centered at cu max ^2Ep 
and is given at sufficiently high temperatures by the 
formula 40 ' - — 



oi(w,j8) 



Biah(§q;/?) (3u> 2 
^wp 16£a 



(4) 
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FIG. 7: Comparison between the fitting curves of the exper- 
imental optical conductivity data at 1.5 GPa using the LP 
model [Eq. ([2]), dashed-dotted line], the SP model in the 
limit s>D [Eq. ((U, dashed line], and the SP model in the 
limit s <C D [Eq. ([5|, dotted line]. The inset displays the 
parameter 2Ep as a function of pressure, obtained from the 
fitting using Eq. ([5]). 



the dc conductivity <r(0, (3) — 34 f2 _1 cm -1 and Ea — 
1331 cm -1 (165 meV). The value of the thermal activa- 
tion energy is consistent with earlier studies^ 4 - However, 
the discrepancy between the fit and the measured data 
is rather large, as seen in Fig. [7] 

For the other limit, s<S^D, the polaronic absorption 
band is calculated according to^ 8 . 



1 



j/T 



-<$>(gj-2E p )N(oj-2E f ), (5) 



where (3 — l/(fcgT), and Ea is the thermal activation 
energy, where Ea=Ep/2. From the fitting we obtain 



<Ji{lo) oc 



where the function $(e) is the corresponding current ver- 
tex which can be derived by $(e) = (D 2 — e 2 )/3, and the 
density of states N(e) — -^jjz^J D 2 — e 2 is assumed to be 
semi-elliptical with the half-band width D. Within this 
model the maximum of the absorption band is positioned 
at uj max — 2Ep — D 2 /2Ep, i.e., it is shifted to lower en- 
ergies relative to 2Ep. The experimental data are well 
described by Eq. ([5]) (see Fig. [?])■ With increasing pres- 
sure the fitting parameter 2Ep is decreased, as displayed 
in the inset of Fig. [71 reflecting the pressure- induced shift 
of the MIR band to lower frequencies. 

According to Fig. [7] the SP model in the limit s<C-D 
captures the shape of the MIR band best, since the fit- 
ting with the LP model gives a too small polaron radius. 
A further test for the quasiparticles relevant for the trans- 
port is the temperature dependence of the dc transport 
data above T v . In case the conduction mechanism is due 
to the hopping of small polarons, one expects a thermally 
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activated behavior of the dc conductivity according tc 44 

T*(T)<xexp[-E H /(k B T)], (6) 

where Eh is the hopping energy. Hereby, the disorder 
energy is omitted, which is justified for crystalline bulk 
materials. The inset of Fig. [T] displays the corresponding 
Arrhenius plot. Obviously, the resulting curve does not 
follow a simple linear behavior. By fitting the data in the 
temperature range 250 - 300 K according to Eq. ^ one 
obtains a hopping energy Eh = 28.5 meV, which is a fac- 
tor of « 5 smaller than the activation energy Ea=Ep/2 
obtained from the optical data within the small polaron 
model (Ep= 271 meV at the lowest pressure applied, 
see inset of Fig. [7]). The discrepancy is, however, re- 
duced with decreasing temperature (see inset of Fig. [TJ, 
indicating that the small polaron regime is gradually ap- 
proached. Interestingly, recent hard x-ray photoemission 
experiments^ also observed a gradual variation from the 
large polaron to the small polaron regime with decreasing 
temperature for the temperature range 250 - 330 K. We 
therefore conclude that the polaronic coupling strength 
in magnetite is intermediate and the prevailing character 
is temperature-dependent. 

For this polaron crossover regime characteristic fea- 
tures are expected in the optical conductivity spec- 
trum, namely peaks at frequencies of the order of the 
phonon frequencies. These peaks correspond to elec- 
tronic transitions between different subbands in the po- 
laron excitation spectrum, denoted as polaron interband 
transitions j 38 i 43 Indeed, the optical conductivity spec- 
trum contains a far-infrared band, which shifts to lower 
energy and whose spectral weight is growing with increas- 
ing pressure (see Fig. [5]). We speculate here about the 
polaronic nature of this far-infrared band. The obser- 
vations are similar to those in other polaronic materials, 
like LaTiC>3.4i and /3-Sr 1 / 6 V2 05 ) 29 ' 33 which were also dis- 
cussed within the intermediate electron-phonon coupling 
regime. 

B. Crossover at ~6 GPa 

Earlier pressure studies found an anomaly in the range 
6-7 GPa in various physical quantities. This anomaly was 
interpreted in terms of a crossover with a rearrangement 
of valence electrons among the different types of Fe sites, 
including the scenario of an inverse-to-normal spinel tran- 
sition according to Eq. (fTJ). In the same pressure range 
(«6 GPa) we observe an enhancement of the pressure- 
induced hardening of the lower-frequency phonon mode. 
But overall, the pressure-induced changes in the opti- 



cal conductivity spectrum are minor. This observation 
is compatible with recent proposals of minor pressure- 
induced alterations of the charge distribution in mag- 
netite, without changing the principle transport mecha- 
nism via polaron hopping motion^ 

A similar pressure-induced rearrangement of charges 
between different transition-metal sites was proposed re- 
cently for /3-Nao.33V205, as revealed by anomalies in 
the pressure-induced shifts of infrared- and Raman-active 
phonon modes? 3 ^^ Similar to the findings in magnetite, 
the charge rearrangement occurs without a change in the 
structural units and the crystal structure symmetry. 4 ? 



V. SUMMARY 

In summary, the optical properties of magnetite have 
been studied at room temperature as a function of pres- 
sure, supplemented by ambient-pressure dc transport 
measurements. The optical conductivity spectrum con- 
sists of a Drude term, two sharp phonon modes, a far- 
infrared band at «600 cm' 1 , and a pronounced MIR ab- 
sorption band, which can be ascribed to excitations of po- 
laronic quasiparticles. With increasing pressure both ab- 
sorption bands shift to lower frequencies and the phonon 
modes harden in a linear fashion. We applied several the- 
oretical models for describing the MIR band and found 
that the SP model in the limit s<C-D captures the shape 
of the MIR band best. However, according to the tem- 
perature dependence of the dc transport data the small- 
polaron character of the charge carriers prevails only at 
low temperature, and therefore the polaronic coupling 
strength in magnetite at room temperature should be 
rather classified as intermediate. This is corroborated by 
the occurrence of a far-infrared band in the optical con- 
ductivity spectrum. For the lower-energy phonon mode 
an abrupt increase of the linear pressure coefficient oc- 
curs at around 6 GPa, which can be attributed to minor 
alterations of the charge distribution among the different 
Fe sites. 
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